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f - "j | Using the generalised set of fluid equations that include the 'Coriolis force' along 

£*\] , with the centrifugal and pressure gradient forces, we have reanalysed the class of 

self similar solutions, with the pseudo-Newtonian potential. We find that the class 
of solutions is well behaved for almost the entire parameter space except for a few 
selected combinations of 7 and a for the co-rotating flow. The analysis of the 

^1^ Bernoulli number shows that whereas it remains positive for co-rotating flow for 

/ > 1/3, for the counter-rotating flow it does admit both positive and negative 
values, indicating the possibility of energy transfer in either direction. 
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1 Introduction 
> 

One of the leading problems, pursued intensely over the last three decades, in 
|^^ , connection with high energy cosmic sources is the topic of accretion disks around 

f~*i ' compact objects. There have been several reviewstrn on the topic covering both 

theoretical and observational findings, but yet the subject is alive with several 
interesting and unsolved problems. In the last j-d|jcade the theory of advection 
C^ . dominated accretion flows drew a lot of attention 000 and this approach is now en- 

Q_i' larging its scope to include convection too in the equilibrium solutions&LJ. As is well 

q , known, the class of solution used for these discussion have been of 'Self Similar' na- 

turc, which looks quite promising in explainintseveral observational featuresuUiij. 
c/3 , Though the classical solution of Narayan & YiH has provided good insight into the 

properties of accretion flows, the solution has a restriction in the parameter space 
with the possible breakdown of self similarity for certain combination of the physical 
parameters, 7 (the gas constant) and a (the viscosity parameter). This was found 
rS , recently by Bhatt & PrasannaEH through an analysis of the solution perturbatively, 

using the well known pseudo-Newtonian potentials. They showed that though self 
similarity could be an excellent approximation over a large range of viscosity pa- 
rameters, for certain combination (a — 0.3, 7 = 1.55), the perturbation develops a 
singularity. In fact it is well known that most of the discussions, particularly in the 
context of modelling accretion around black holes and neutron starsu'EJ'tLa, assume 
the value of gas constant 7 to range between 4/3 to 5/3. Since it was shown that 
a singularity could develop for a value in between, it is indeed necessary to find 
solutions without the pathology at least in the parameter regime of interest. 

It is also known that almost all the discussions, particularly in the Newtonian 
regime, consider the equations of motion for the fluid in a background geometry 
where the effects of rotation of the central object is not considered. In the case of 
black hole accretion, a correct model should be fully general relativistic, which if 
considered on a Kerr background geometry would automatically take into effect the 
rotation of the black hole through 'inertial frame dragging' effect. If one wants to 
create a similar scenario in a purely Newtonian discussion the only way would be 
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to bring into focus the 'Coriolis force', wherein one writes the equation of motion 
in a rotating frame. In fact the concept of reintroducing the language of "incrtial 
forces" in general relativity, though not very popular, has beenJjcuindJxi be useful 
in explaining some already known results more satisfactorily EjiijOa. Particu- 
larly the fact that the "Centrifugal force" reverses sign very close to the compact 
object, could be a plausible reason for some of the known effects. As was shown 
by Abramowicz et al.Ej it is indeed possible to introduce the concept of inertial 
forces in a covariant formalism that is valid in general spacetimes, without any 
particular symmetry requirement. Of these forces, the "Coriolis" force is generally 
assumed to be of less relevance in most of the Newtonian physics. However, in 
general relativity the term which gives rise to the Coriolis rt¥pe field is the well 
known "Lense-Thirring" effect of dragging of inertial frameset, which couples the 
angular momentum of the gravitating source with the particle dynamics, through 
the spacetime curvature. It is indeed interesting to note that a particle in circular 
orbit around a Kerr black hole at a radius wherein the centrifugal force is zero, 
has to change its angular velocity depending upon the Kerr parameter 'a' to be in 
equilibrium. In fact while a co-rotating particle has to decrease its angular velocity 
slightly with increase in 'a'. tb p counter-rotating ones have to increase it substan- 
tially to keep in equilibriumliilEl. Though this is purely a general relativistic effect, 
which happens very close to black hole, one wonders whether in principle there 
could be effects of the rotation of the central gravitating source on the particles or 
fluids at a distance, even in a purely Newtonian approach. It is quite well known 
that while discussing disks in binary systems and planetary rings one does consider 
the equations of motioiL m a rotating coordinate system, that includes the Coriolis 
and centrifugal terms Etlj. 

With this background we now look for possible influence of Coriolis type terms 
on the fluid flow in accretion disks around compact rotating objects, particularly for 
the type wherein the fluid disk is in contact with the central body. The present aim 
of this discussion is to find out whether the introduction of the rotational effects 
through 'Coriolis' type terms in the equations of motion would alter the situation 
regarding the occurrence of singularity in the perturbation solution. Further, as 
will be shown, this treatment also suggests a way to distinguish between the co and 
counter rotating accretion flow onto compact objects. 



2 Formalism 

As is usual in accretion disk theory we consider the height integrated set of equa- 
tions describing the steady state, axisymmetric fluid distribution that allows one 
to describe all physical quantities as functions of the cylindrical radius R only. 
Mass conservation equation yields 

- inRHpV = M = constant , (1) 

wherein H = ■JW/2(c s /Q,k), is the vertical half thickness with c s and VLk being the 
isothermal sound speed and the Keplerian angular velocity respectively. 

The equation of motion for the fluid in the gravitational field of a slowly rotating 
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compact object is given by 

^ + 2l^xf + ^x^x]? = -^ + l t „ + I /V 2 f (2) 

dt p 

wherein, ~uf is the intrinsic angular velocity of the central body, V the 3-velocity 
of the fluid element, ~F g the gravitational acceleration and v is the coefficient of 
kinematic viscosity — ac 2 /Qk n - As we are using the height integrated equations we 
get from (2) the radial and angular momentum conservation equations 

y^L - m 2 - 2^R - w 2 R = -m 2 --^- (3) 

dR kn pdR w 



and 



v m^ = ^mm{» pFi ' H m)-^ RV < 4 » 

with V and £1 being the radial velocity and angular frequency respectively of the 
fluid element. 

We shall use the quasi-Newtonian potential 

4> = GM/{R -R g ), R g = 2GM/c 2 (5) 

in terms of which the Keplerian angular velocity is 

il 2 n = GM/R(R-R g ) 2 . (6) 

As mentioned earlier we shall now couple w and fi through a simple relation 
lo = aQ, a being a constant. Using the equation of state p = pc 2 we can rewrite (3) 
and (4) as 

"£-*»" + n ^ + *-fg + 7ffi- (7) 

with n = (1 + a) 2 and 

ffi-sS*[ n * + '/ anw H (8) 

where the constant j has the usual interpretation of being the angular momentum 
per unit accreted mass. However for self similar solutions one assumes j negligible 
in comparison with the other terms. 

Apart from these hydrodynamical equations one also considers the thermody- 
namical balance between the local viscous heating and radiative cooling that gives 
rise to advection and is obtained through the energy equationB 

pV d^ 2 dp_ fapc 2 R 2 fdn\ 2 



-y-ldR s s dR fl kn \dR 
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3 Equilibrium Solutions 



The original set of self similar solution were obtained by setting -^ = in the po- 
tential, thus taking Qk n = Qk = (GM/R 3 ) 1 / 2 and assuming the radial dependence 
of the other parameters to be 



V = V 



R a \ 1/2 (B.^ 1/2 



= n Q 



3/2 ,„ x 3/2 



Assuming the same dependences and substituting these in Eqns. (1), (7), (8) 
(9), one can get the set of algebraic equations 



V 2 - 1 + 2nQ 2 Rl + 54 = 



2o 



2 C fo = -^o(l + 4a 



(37-5)V = ^=(7-l)/a<W 



po = -M/(4irV5c s0 V R 



(10) 



Solving these, one can get 



V = -Ag/(3y/2a), 

c 2 sQ =Ag(l + 4a) /9a 2 

n R g = (e 1 'Agf' 2 /3a, 



Po 



63 



-(Agf(l + 4a) 



1-1/2 



(11) 



with 



(I 



A = 2nd 



-1 



(7- 
5(1 



7) 



1)/ 
-4a) 



1 



18a 2 
~A^ 



(12) 



Figures 1 and 2 show the behaviour of the various physical parameters for the cases 
a > and a < respectively. The behaviours are shown as a function of R for 
fixed values of a and 7. The general trend of the solutions are similar for all values 
of a (to be expected) except for the fact that with a^O the relative values tend to 
decrease with higher a values for the fluid velocities ( Fig. lab) in the co-rotating 
case, but increase for the counter-rotating flow (Fig. 2ab). Whereas the sound 
speed shows no variation for the co-rotating case (Fig. lc) but very small change in 
the counter-rotating case (Fig. 2c), the density shows an increase with the values 
of a for the co-rotating flow (Fig. Id) and decrease in the counter-rotating cases 
(Fig. 2d), the change being pronounced only close to the compact object. 
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Fig. 1 
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Fig. 1 : Plots of radial velocity V, angular velocity Q, sound speed c s and density p as 
a function of R for different values of the rotation parameter 'a' for fixed a = 0.3 and 
7= 1.5; o = ( ),o = 0.5( ) and o = 1.5 (- ----). 



Perturbation 



Following the approach of Bhatt & PrasannaO we now study the effects of non- 
Newtonian potential on the self similar solutions considered above, through pertur- 
bations 

V — * Po + 8p, V — ► Vq + SV, fi — ► Qq + Sfl, c s — > c s0 + Sc s 

in the equations of motion and mass and energy conservation equations. We then 
simplify by retaining the terms linear in the perturbations, after using the back- 
ground solution (11) and obtain the equations governing the perturbations. Per- 



turbation in the Keplerian velocity to the first order in -£ 

sn 2 kn « 2^n 2 K , sn kn = ^n K 

1 /2 

where CIk is taken as (GM/R 3 ) 

Thus the set of equations governing the perturbations 



is given by 
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Fig. 2 : Plots similar as Fig. 1 but for a — (- 
-0.2 ( ). 



= —0.1 ( ) and a = 



SV, <5f2, £c s and Sp are: [] 



<5# Sp SV SH _ Sc s R g 



.dV 



Sp dp Id 



V dB 6V + 5V ~ rB = 2nRQm ~ Rm kn + -31^- -Tn S P 



dR "" m ' p 2 di? pdiT 

<5p = Cg(5p + 2pc s Sc s , (13) 

n + % f siRdRj {n kn sv + vsn kn )+vn kn un + J /" am^r J , 



2 d dfl 

ac s -—5\l+2ac s -—Sc s = 
art art 

/^ rfcg 
7-ldi? 

a There are a couple of minor typographical errors in liil, which can be corrected by taking these 
equation with a=0. 



Sp SV n fdc 2 s \ d . . , 



s dR 



Sc s SV (d P \ l d 
2 ^7 + - + {dR dR Sp 
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faR 2 c 2 s p (dQ. 



<1R 



Cti 



sp 2^ _r 1 + 2 (dn\~ 1 ± sn 

p c s R \dR dR 



Here again as in the equilibrium case, for self similar solutions, the assumption 
j << S(flR 2 + 2a j VtRdR) is maintained. Thus again using self similarity of the 
perturbed solutions through 



SV = Vi,i{ -R 



3/2 



, Sc s = c s0 i [ -^ 



3/2 



■ «> = <> lU (^ 



5/2 



Rn 



^ = Poi ( ^ ) (14) 



and substituting these in (13) one gets the set of algebraic equations: 



5/2 



9o 



^(^-(1 



*< 



2 Vo 



Po 






CsOl 
C S 



1, 



2ne'Ag fi i 6^(1 + 4a) c s0 i 
9a 2 fin + 



3(1- 4a) 

2(1 + 4a) 



9a 2 

"^7 



- 1 



CsO 



, c s01 



Ag(l + 4a) 
9a 2 



(15) 

= 0, 

(16) 

(17) 



(7- 



»t 



10 



|-7)^-2(7 + l)^=7-5 7 . (18) 

3 / Ho c s0 



Since it is a consistent set of inhomogeneous equations, one can solve for the per- 
turbations and thus get the set of solutions incorporating the quasi Newtonian 
potential, to be 

3/2 



P P = Pa 



V P = V 



CsO 



Rg 

R 

Rg_ 

R 

Rg 
R 

fi» = On I —— 



1 



1/2 



1/2 



3/2 



Poi 

Po 


Rg 

R 


v Q1 

V 


Rg 

R 


CsOl 


Rg 



(19) 



Cso R 

^01 Rg 

no~~R 



i + ^^f 



Figures (3)- (5) present these solutions for the physical quantities as a function 
of R with two of the three parameters a, 7, a fixed, while the other is varied. While 
the changes as a function of a remains almost the same as in the non-perturbed 
case (Fig. 3), there appear some changes as a function of the viscosity (a) and gas 
parameter (7). As may be noted in Fig. 4a the radial velocity shows a jump in the 
profile as 7 changes from 1.1 to 1.6, while the other three do not show such a jump 
(Fig. 4b,c,d). On the other hand as is depicted in Fig. 5, their behaviour as a 
function of a for fixed 7 = 1.5 and a = 0.5 is quite interesting. While V p increases 
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Fig. 3 
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Fig. 3 : Same as Fig. 1 but with the addition of perturbative solutions arising from quasi- 
Newtonian potential; a = ( ), a — 0.5 ( ) and a = 1.5 (— — — — — ). 



with a (Fig. 5a) and p v decreases (Fig. 5d), £l p and c p show no change at all with 
varying a (Fig. 5b, c). 

In order to understand these features, we have presented the plots of the ratios 
of perturbed to the unperturbed quantities (Voi/Vo, Ooi/fio> c s qi/c s o, Poi/po) in 
Figs. (6) and (7) one as a function of 7 and the other as a function of a for the 
cases a = 0, 0.5, 1.5. As is evident from Fig. 6 the singularity in the perturbations 
has shifted to lower values of 7 as a changes from to 1.5. For a = 0, the location 
of the 

singularity in the 7 space is at 7 = 1.55, a result which was obtained earlier by 
Bhatt & PrasannaO. Figure 7 shows the perturbation in the a-space, whereas for 
a = the perturbations grows as the flow gets nearer to the compact object, (a 
result known earlier), with the inclusion of Coriolis terms a ^ 0, the perturbations 
are well confined and further the ratio decreases as the flow approaches the central 
gravitating source. Fig. 8 depicts the plots of the ratios of the perturbed to 
unperturbed physical quantities for a = 0, —0.1, —0.2 as a function of 7 for fixed 
a = 0.3. It can be clearly seen that there appears no singular behaviour for a ^ 
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Fig. 4 
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Fig. 4 : Same as Fig. 3 for fixed a — 0.3 and a — 0.5, but with different 7s; 7 = 1.1 

( ), 7 = 1-4 ( ) and 7= 1.6 ( ). 

unlike the case a = 0. 



5 Discussion 

The introduction of the Coriolis force and a direct coupling of the two angular 
velocities seem to bring in some features which could be helpful for understanding 
the parameter space available for model building. At the outset one can see that the 
reality of the solutions, immediately requires (l+4o) > 0, meaning that the counter- 
rotating flows have to have a minimum speed, at least greater than four times the 
angular velocity of the central source, to be in equilibrium. This may not be very 
surprising, as the influence of the central rotation of 'frame dragging' would imply 
that any counter-rotating body has to overcome this barrier. Though in principle 
'frame dragging' is purely a general relativistic effect, it appears as though the 
Coriolis term in Newtonian physics mimics a similar effect. An important outcome 
of the analysis is that the rotational effect seems to stabilise the flows for the usually 
considered gas and viscosity parameters 7 and a respectively as depicted by the 
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Fig. 5 : Same as Fig. 3 for fixed 7 = 1.5 and a — 0.5, but varying a; a = 0.01 ( ), 

a = 0.1( )anda = 0.3(- ----). 



perturbations of the self similar solution; (4/3 < 7 < 5/3, < a < 0.3). 
Equation (7) shows that the solution satisfies the relation 

Yl 

2 



— + nn 2 R 2 - Vl\R 2 + -cl = 0. 



(20) 



Using this in the computation of the normalised Bernoulli parameter '6' as defined 
by Narayan & Yi H 



1 (W &# 2 2 

b -v 2 {— + — - nKR T ~ 



7 c? 



one finds after using Eqns. (11) and (12) 

Age' 



6=-^[l-2n + 3/(l + 4a)] 

19a z 

(l-2n + 3/(l + 4a)) 
hn + ^ (1 + 4a)) 



(21) 

(22) 
(23) 
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Fig. 6 
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Fig. 6 : Plots of the ratios of perturbed to the unperturbed parameters as a function 
7, for fixed a = 0.3 and different values of a; a = ( ), a = 0.5 ( ) and 

a = 1.5 (- - - - -). 



Figures (9) and (10) show the behaviour of '6' for different ranges of parameters a, 
7, /. For 1 < 7 < |, b is positive for / > 1/3, for a < —0.1 and as a approaches 
—0.2, b tends to become negative even for values of / > 1/3. As a tends to the 
value —0.25, b is negative for all values of / and 7. Thus while the co-rotating 
fluid can have energy transfer only outwards for / > 1/3, the counter-rotating fluid 
can have it in either direction depending upon its angular velocity as compared 
to that of the central accreting source. For advection dominated flows (/ = 1) b 
is mostly positive for a > —0.2. However b changes sign for co-rotating flows at 
a = 2 + V5 and for counter-rotating flows at a = 2 — v5- Hence in principle if the 
energy transfer inwards has to be effective, the co-rotating flow has to have very 
low angular velocity Q < w/(2 + v / 5), whereas the counter-rotating flow has to have 
very large angular velocity f2 > lo/ (2 — \/5)- 

In conclusion it may be seen that introducing rotational effects into a Newtonian 
description of the accretion flow, through the Coriolis term, has indeed enlarged 
the parameter space of the self similar solutions, making them viable for model 
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Fig. 7 
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Fig. 7 : Plots same as in Fig. 6 for fixed 7 = 1.5 as a function of a and different values 
of a; a = ( ), a — 0.5 ( ) and a — 1.5 (— — — — — ). 

building scenarios, particularly treating both co and counter rotating flows. A sim- 
ple coupling of the angular velocities has distinctly shown the possibility of energy 
transfer both inwards and outwards depending on the effective angular velocity of 
the fluid flow. More detailed studies, like looking for a global solution and inclusion 
of other effects like convection, may perhaps reveal interesting features which could 
influence the 'energy budget' of accretion dynamics. 
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